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INTRODUCTION 


The high temperature environment of cruise flight near Mach 3 makes vane-type 
flow direction sensors generally inadequate for use on the YF-12 aircraft. As a 
result, angle of attack and angle of sideslip are measured by a fixed four-port 
pressure-sensing hemispherical head mounted on the nose boom . The presence of 
time lag errors in pitot-static instrumentation is documented in references 1 to 3 , in 
which it is noted that the lag can be as great as several seconds at high altitudes 
and airspeeds. Less attention has been paid to lag in angle-of-attack and angle-of- 
sideslip pressure sensors . The lag of these sensors is considerably less than that 
in pitot-static systems; however, a small lag in the measurement of angle of attack 
and angle of sideslip can be critical when these variables are used in stability 
derivative determination or as inputs to a control system . For example , in the 
YF-12 flight program, a lag in the sideslip angle sensed by the inlet control system 
caused the inlet forward bypass door to drive the vehicle unstable (ref. 4) when 
the stability augmentation system (SAS) was switched off. 

Theoretically lag has been determined by computing estimates based on viscous 
lag theory (as in ref. 2) and experimentally by applying known pressure variations 
to the instrument pressure ports in laboratory ground tests . This report presents 
three other methods of determining lag from dynamic stability and control maneu- 
ver data: time vector, Newton-Raphson, and least squares. Lags of the YF-12A 
angle-of-attack and angle-of- sideslip sensor systems are documented throughout 
the flight envelope , and the results are compared with those obtained by using 
viscous lag theory . 


SYMBOLS 


Physical quantities in this report are given in the International System of 
Units (SI) and parenthetically in U.S. Customary Units. The measurements were 
taken in U.S. Customary Units. Factors relating the two systems are presented in 
reference 5. 
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normal acceleration at center of gravity , g units 

transverse acceleration at center of gravity , g units 
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acceleration due to gravity, meters/second (feet/second ) 
pressure altitude , meters (feet) 
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pressure, newtons/meter (pounds/foot ) 


roll, pitch, and yaw rate, respectively, radians /second (unless 
otherwise noted) 

2 2 

dynamic pressure, newtons /meter (pounds /foot ) 

Laplace operator 

temperature , degrees Kelvin (Rankine) 
time, seconds 

velocity, meters/second (feet/second) 

distance from center of gravity to angle-of-attack and angle-of- 
sideslip sensor , respectively , positive forward , meters (feet) 
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Subscripts: 

c 


angle ol sideslip at center of gravity , radians (unless otherwise 
noted) 

ratio of specific heats ,1.4 for air 
incremental change 

differential pressure measured by the sensor in the a and p 

2 2 

plane, respectively , newtons /meter (pounds/foot ) 

aileron , elevator , and rudder deflection , respectively , radians 
(unless otherwise noted) , positive when aileron deflection 
produces right roll , elevator trailing edge down , trailing edge 
of rudder left 

damping ratio 

Euler angle of pitch and roll , respectively , radians (unless other- 
wise noted) 

2 

coefficient of viscosity , newtons-second/meter 

2 

(pounds-second/foot ) 
standard deviation 
time constant , seconds 
damping angle, sin *(£), degrees 

phase angle of quantity j relative to quantity k , degrees 

damped and undamped natural frequency of the airplane in 
oscillation, respectively, radians/second 


calculated 


f flight 

i indicated 


o initial 
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stagnation condition 


free- stream condition 

angle of the a and p ports on the hemispherical head with respect 
to the centerline of the probe , degrees 

port angle from centerline of hemispherical head 

A dot over a quantity denotes the time derivative of that quantity . 


airplane AND INSTRUMENTATION 

The YF-12A airplane (ref. 4) is a twin engined, delta-winged interceptor 
designed for sustained cruise above Mach 3 and altitudes greater than 24,400 meters 
(80,000 feet) . As noted previously, the high temperature environment of Mach 3 
cruise flight necessitated the use of a fixed four-port pressure sensing hemispherical 
head to obtain angle-of- attack and angle-of-sideslip measurements on the airplane. 
Reference 6 describes the a-p sensor and its calibration in detail. The sensor , 
attached as a dogleg to the nose boom probe , is shown on the airplane in figure 1 . 


Test airplane with nose boom probe 

E-23131 


Figure 1 

attached 



A closeup of the assembly is shown in figure 2, and a dimensioned three-view 
drawing in figure 3. 


boom probe showing the a ~^ E ^g°g 


Figure 2. Closeup of nose 
assembly. 


Standoff 

distance 

J2.50I 


Port offset, 3.4* 
from centerline 


Figure 3. Nose boom probe dimensions fin centimeters finches)) with 
details of sensor head configuration. 



Basically , the a-p sensor provides differential pressures , AP and AP Q , 
as shown in the following sketch: ” 



The sensor is connected to pressure transducers in the nose cone of the air- 
plane . Geometric characteristics of the system are shown in the following schematic 
drawing and tabulation: 


First segment Second segment 



Total volume of 
ail transducers 

Sensor system 

a 

p 

Port - 



Depth ♦ cm (in . ) 

0.203 (0.080) 

0.203 (0.080) 

Inside diameter , cm (in.) 

0.178 (0.070) 

0.178 (0.070) 

Tubing - 



First segment: 



Length, cm (in .) 

43.2 (17) 

43.2 (17) 

Inside diameter , cm (in.) 

0.365 (0.144) 

0,365 (0.144) 

Second segment: 



Length, cm (in.) 

640 (252) 

640 (252) 

Inside diameter , cm (in.) 

0.493 (0.194) 

0.493 (0.194) 

Transducer volume - 



a or p pulse code modulation (PCM) instrumen- 



3 3 

tation transducer . cm (in ) 

16.4 (1) 

16.4 (1) 

Sum of remaining transducers used for onboard 



3 3 

systems .cm (in ) 

239.2 (14.6) 

157.3 (9.6) 

Total 

255,6 (15.6) 

173,7 (10.6) 


Wind-tunnel calibration curves of (AP Q /q)/a and (APp/q)/p as a function of 

Mach number for this probe were obtained from reference 6 and are shown in 
figure 4. These curves were applied to the AP Q /q and APp/q flight data to obtain 
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values of a and p as a function of time. For the maneuver analyses, dynamic 
pressure was assumed to be constant and was calculated from the initial Mach 
number and altitude conditions . 


(AP I qUa, (APJ qUP, 
a p 

per deg 



Figure 4. Calibration curves for conversion of differential 
pressures to angle-of-attack and angle-of-sideslip values 
as a function of Mach number. (Adapted from ref. 6.) 


Basic stability and control instrumentation provided angular rates , normal and 
lateral accelerations , and pitch and roll attitudes . 

Pertinent characteristics of the PCM system instrumentation used for the time 
lag analysis are listed in the following tabulation: 


Parameter 

System 

range 

System 

resolution 

to , rad /sec 
n 

Pressure, P 

±20,261 N/m 2 

79 N/m 2 

3600 


(±424 lb/ft 2 ) 

(1 .65 lb/ft 2 ) 


Roll rate, p 

±20 deg/sec 

0.078 deg/sec 


Pitch rate, q 

±10 deg/sec 

0.039 deg/sec 

■ 

Yaw rate, r 

±10 deg/ sec 

0 . 039 deg/ sec 


Euler angle of roll , <p 

±60° 

0.234° 


Euler angle of pitch , 0 

-15°, 25° 

0.078° 


Normal acceleration , a 

n 

-lg. 3(? 

0 . 0078g 

500 

Transverse acceleration , a^ 

±0.5* 

0 . 0020g 

500 


Data were recorded on magnetic tape by means of the PCM system and were analyzed 
at a rate of 40 points per second . 


TEST MANEUVERS 


The data used in this analysis were obtained in a previous study for which 
standard stability and control maneuvers were performed . Longitudinal data were 
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from elevator pulses, and lateral-directional data were from rudder and aileron 
doublets . The stability augmentation system was not used during the maneuvers so 
that damping would be as low as possible. 


DATA ANALYSIS 


Three methods of analysis were used to determine the lag of the a and p 
systems: time vector, Newton-Raphson, and least squares. For the purpose of 
this report, the o and p sensor systems are considered to include the probe, the 
tubing, the transducers, and the data acquisition system. 

The time-vector method is suitable for use in analyzing free-oscillation maneu- 
vers with damping ratios less than 0.3. Time histories calculated with the Newton- 
Raphson technique can be used to define lags for maneuvers that are oscillatory but 
not necessarily completely free of control inputs . The least-squares method , the 
most general of the three , is not restricted by high damping ratios or control inputs . 
With all three methods , the lag was determined by comparing flight data with calcu- 
lated time histories . 


Time-Vector Method 


The time-vector method (ref. 7) establishes amplitude ratio and phase angle 
relationships between the different parameters in the linearized equations of motion. 
By using the appropriate vector format equations , shown in appendix A , phase 
relationships of a and p to other flight variables can be established without using 
a and p. 

For a longitudinal maneuver (fig. 5) with pitch rate, q, and normal accelera- 
tion, a n , available, the phase angle, ^ , can be calculated. The lag in the 

indicated a shows up as the discrepancy between this computed angle and the 
flight data, 


For a lateral-directional maneuver , a similar procedure using roll rate , p , yaw 
rate, r , and transverse acceleration , a^ , produces a calculated *** (p./ r ) • The 

difference in phasing between the calculated data and the flight data, 4 > (n , can 
be attributed to the lag in the indicated p . 


(Pj/ry 


The requirements of low damping and free oscillation limit the number of maneu- 
vers suitable for time-vector analysis , and the manual nature of the method makes 
analycis of low-quality data (for example , noise-contaminated data) susceptible to 
analyst-induced biases . 
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Figure 5. Typical time history of a longitudinal maneuver suitable 
for time-vector analysis. M = 2.78; h = 20,730 meters (68,000 feet) . 

Newton-Raphson Derivative Matching 

The basic principle of the Newton-Raphson derivative determination technique 
(ref. 8) is the minimization of the error in matches of flight and computed time 
histories. Use of this method results in a complete set of stability and control 
derivatives , as well as a calculated time history of all parameters being matched . 
Lags in any of the parameters , if not taken into account , will degrade the quality of 
the match by interfering with the overall phasing . Appendix B presents the equa- 
tions used in the Newton-Raphson technique. 

In the lateral-directional mode , matching p is not required to obtain a solution , 
which makes this method suitable for lag determination. A Newton-Raphson match 
of time histories of a supersonic lateral-directional maneuver is shown in figure 6 . 
The phase angle between the calculated and the flight angle-of-sideslip traces gives 
an approximation of the lag magnitude . A more accurate method of determining the 
lag is to compare the phase angle of the flight data, / r ) > with tne phase angle, 

0 . . , determined from the calculated time history and assign the difference to 

( (V r) c 

lag in angle of attack. In this way, by comparing the flight and calculated values 
with a primary flight parameter , any possible overall phasing discrepancy of the 
match as a whole- is taken into account. 
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Figure 6. Newton-Raphson match of time histories of a lateral-directional 
maneuver . M =2.76; h = 20,910 meters (68,600 feet) . 

Our experience has shown that the Newton-Raphson method is not well suited to 
determining lags in data from longitudinal maneuvers. In addition to normal accel- 
eration, pitch rate, and pitch attitude, the angle-of-attack trace must be used to 
arrive at a successful Newton-Raphson match. Using the uncorrected flight a will 
result in a match , but the derivatives are thereby biased , and the quality of the 
match , as measured by the fit error of the parameter time histories , is not the best 
obtainable. In this instance the a trace must be adjusted for lag to obtain a reliable 
analysis. 





Least-Squares Fit Error Minimization 

The least-squares technique is a general technique in that time histories of all 
types of maneuvers are suitable for lag determination . Time histories of a or p 
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for a given maneuver were calculated from accelerations and angular rates using 
the equations in appendix C. Figure 7 presents flight and calculated time histories 
of two longitudinal maneuvers and one lateral -directional maneuver . Figure 7 (a) 
presents the a traces for the longitudinal maneuver shown in figure 5 , and 
figure 7 (b) presents p for the supersonic lateral-directional maneuver shown in 
! ,mre 6. The highly damped a response typical of subsonic longitudinal dynamics 
i. 'hown in figure 7(c) . 



t, sec 

(a) Longitudinal maneuver. M = 2.78; 
h = 20, 730 meters ( 68, 000 feet) . 



t, sec 


(b) Lateral -directional maneuver . M = 2.76; 
h = 20,910 meters (68,600 feet) . 



t, sec 

(c) Longitudinal maneuver. M = 0.80; 
h = 6, 460 meters ( 21 , 200 feet) . 

Figure 7. Flight and calculated time histories of angle of 
attack and angle of sideslip for three maneuvers. 
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Plotting a calculated time history against a flight a or p trace at some chosen 
sample rate over the same time segment would produce a straight line if there were 
no time lag. The presence of a lag is shown as a phase shift between the two time 
histories and introduces a loop into the plot . Cross-plotting the data of figure 7 
results in the ellipse type of plot shown in the upper half of figures 8(a) to 8(c) . 
The standard deviation, o, from a least-squares straight-line fit through the points 
of the loop provides a quantitative indication of the time lag. 
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- 2.00 - 1.60 - 1.20 -.80 -.40 0 . 40 . 80 1.20 1.60 2.00 

a c , deg 

(a) Supersonic longitudinal maneuver shown in figure 7(a). 
Standard deviation minimized for At = 0.40 second. 

Figure 8. Basic and time-shifted cross plots of flight and calculated angle 
of attack and angle of sideslip for three maneuvers shown in figure 7. 
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- 2.00 - 1.60 - 1.20 -.80 -.40 0 . 40 . 80 1.20 1.60 2.00 


P £ , deg 

(b) Supersonic lateral-directional maneuver shown in figure 7(b) 
Standard deviation minimized for At = 0.35 second. 

Figure 8. Continued. 









(c) Subsonic longitudinal maneuver shown in figure 7(c). 

Standard deviation minimized for At = 0.15 second. 

Figure 8. Concluded. 

As the flight time history of a or p is shifted backward in time with respect 
to the calculated a or ? , the lag is reduced and the size of the loop decreases , as 
shown in the lower plots of figures 8(a) to 8(c) . The time shift, At, that minimizes 
o defines the magnitude of the lag within the resolution of the chosen sample rate 
of 40 points per second . The trend of o with At for a i three maneuvers is shown 
in figure 9, in which the corresponding minimums of standard deviation are well 
defined . 
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o a. M = 2.78, h = 20.730 m (68,000 ft) 



Figure 9. Variation of standard deviation as a function 

of time shift for the three sample maneuvers . I 

i 


Even though the least- squares method is a general technique, the best results j 

are obtained from analyses of highly dynamic maneuvers and maneuvers with low j 

effective damping . j 

I 

ACCURACY ! 


The application of the Newton-Raphson and the least-squares techniques pro- 
duces solutions with minimum errors , as stated by the definition of their mathemat- 
ical development . The minimum error is in turn dependent on the accuracy of the 
mathematical model used and the quality of the data . An exact assessment of the 
effect of model and data errors on the solution is complex and beyond the scope of 
this report . However , it is estimated that the error in the measured time lag is 
approximately ±0.05 second. This estimate is based on the high accuracy and fre- 
quency response of the data system instrumentation, experience with stability and 
control analysis using similar data systems , and a review of the lag results . It is 
estimated that the combined effects of the PCM system , parameter calibration , and 
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various instrument characteristics produced less than 1 -percent full-range error 
on the parameters used . In perturbation equations , such as those in this report , 
the constant or bias errors canceled out , and the parameter errors were less than 
1 percent . Figure 9 illustrates , qualitatively , the uniqueness of the least-squares 
solution . 

The influence of aircraft flexibility on the results was assessed and found to be 
negligible . 

The time-vector approach requires a visual fairing of the data to obtain ampli- 
tude ratios and phase angles . Therefore , the accuracy of the results is similar to 
that of the Newton-Raphson and least-squares methods. 


RESULTS AND DISCUSSION 


Flight-Determined Lag 
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The a and p sensor system configurations were essentially the same; there- 
fore , no marked difference between lag in a and p was anticipated . Lag results 
from all three methods are 
presented as a function of 
Mach number and knots 
equivalent airspeed in 
figures 10 to 12. 

Time-vector-determined 
lag results are presented in 
figure 10 for the angle-of- 
attack and the angle-of- 
sideslip systems. An appar- 
ent increase in lag with 
increasing Mach number 
is shown; however , the 
increase ir. Mach number 
also generally represents 
increasing altitude . The 
data are not sufficient nor 
definitive enough for Mach 
number and altitude effects 
to be separated . For the 
a system , only supersonic 
data are shown because 
high damping precluded 
analysis of data from sub- 
sonic maneuvers. The 
scatter of the data is 
nominal , inasmuch as 
the sample rate interval 
was 0.025 second. 
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Figure 10. Lag determined by time-vector method 
as a function of Mach number and knots equivalent 
airspeed. 
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Newton-Raphson-determined A ag for the p system is shown in figure 11. It 
was not possible to determine angle-of-attack lag with this technique because angle- 
of-attack matching is required to obtain consistent Newton-Raphson results. The 
lag trend with Mach number is similar to that obtained with the time-vector method , 
and the variation is approximately equal to the PCM sample rate interval. 
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Figure 11. Lag determined from Newton-Raphson 
matches as a function of Mach number and knots 
equivalent airspeed. 


Least-squares-determined lag results are shown in figure 12 for the angle-of- 
attack and the angle-of-sidesli.> systems. The trends are the same as those deter- 
mined by the time-vector and Newton-Raphson techniques. There are many more 
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Figure 12. Lag determined by least-squares technique as a 
function of Mach number and knots equivalent airspeed . 
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data points than for the other techniques , because the general nature of the least- 
squares technique permits analysis of all types of maneuvers. The automatic nature 
of the technique eliminated analyst-induced errors and reduced data scatter . 

Although all three methods produced similar results , the least-squares tech- 
nique was superior because of its general applicability and the improved data con- 
sistency , thus only the results shown in figure 12 are used in the following discus- 
sion. 


Conversion of Lag to Time Constants 

Most literature on pressure lag assumes a first-order lag model with a Laplace 

transfer function of the type -r -~ — , in which t is the time constant . For meaning- 

ful comparison of the At lag values with existing predictions , the values must first 
be converted to t, since At is frequency -dependent in this model. With the data 
available , this was done by using the following equation , developed in reference 4: 

-1 

... . 05 d 

^n + tan p7.3(At)o> d J 

The frequency and damping parameters are those of the oscillation and were calcu- 
lated from the stability and control derivative results. 

Inasmuch as altitude instead of Mach number has traditionally been considered 
the major correlation parameter for lag (refs . 3 and 9) , a plot of the least-squares- 
determined time constants for a and p as a function of altitude is shown in 
figure 13. The Mach number for each data point is indicated. Because of the 
limited flight conditions under which the stability and control maneuvers were 
performed, Mach number and altitude are not entirely independent in these data. 

An apparent trend of increasing lag with increasing altitude is shown , but the 
correlation between Mach number and altitude makes it impossible to separate their 
effects on the lag . 


Theoretical Time Constants 

Many methods of calculating theoretical time constants of pressure sensing 
systems have been developed. Three methods in particular were considered in 
this study in attempting to correlate flight results with theory. In the first method, 
which could be regarded as the classical technique , the system must be repre- 
sented by only one tubing segment and one transducer . The second method 
considered was developed by Lamb (ref. 10) and is basically an extension of the 
previous method in that a variety of geometrical configurations can be considered . 
The third method was developed by Iberall (ref. 11) and is based on oscillatory 
pressure theory . 
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Figure 13. Least-squares-determined time constant plotted 
against altitude. Mach number is indicated for each data point. 
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The results from all three methods were similar for the system under consider- 
ation; therefore, only the classical theory results are presented. Spot check 
calculations using the other two methods were made , and the manner in which the 
results would differ from those of the classical theory is discussed . 

In the classical theory , the time constant , t , is proportional to the ratio of 
p/P , in which p is a function of the system temperature and P is the pressure , 
P 450 , sensed at a port on the hemispherical head. Appendix D presents the method 

of calculating the classical theoretical time constant using only the dominant segment 
of tubing and the sum of the transducer volumes. Figure 14 compares the theoretical 
time constant of the a system with flight data. (The theoretical p system lag is 
22 . 5 percent smaller than the a system lag.) There is essentially no agreement, 
because the flight-measured lags are approximately four times as great as the 
theoretical lags . 


0 10 20 30 40 50 60 70 80 x 10 3 

1 1 1 1 1 1 1 1 1 



Figure 14. Comparison of flight time constants for the 
a and p systems with theoretical time constants of 
the j system. 


Theoretical calculations based on the method developed by Lamb indicated that | 

tie port and the first 43.2-centimeter (17-inch) segment of tubing would increase ! 

the lag shown in figure 14 by 28 percent , which is still much lower than that for the i 

flight data. Calculation '> made using the oscillatory pressure theory developed by I 

Iberall resulted in lag values similar to those obtained by using Lamb's technique. 1 

However , it cannot be concluded that the classical theory is incorrect . The a and 1 

P system lags were calibrated in the laboratory , and the results agreed with those | 

of the classical theory . In the laboratory , known pressure variations (primarily 1 
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ramps) were applied directly to the ports, whereas in flight the flow field around 
the hemispherical head determined the port pressure. The results in this report 
are based on dynamic maneuvers; hence a and p changed with time, which in 
turn caused a mass flow in and out of opposing hemispherical head ports . A primary 
contributor to the lag discrepancy could be the mass flow interference around the 
port, which in turn would affect the hemispherical head pressure field. Reference 12 
discusses this flow interaction phenomenon as it relates to static-pressure measure- 
ments in climbs and descents for a typical pitot-static nose boom installation . Other 
factors not accounted for in the simple viscous lag theory are the dynamic flow 
effects of compressibility, mixed flow, and heat transfer (ref. 13) . 


CONCLUDING REMARKS 


Flight data from the YF-12A airplane were analyzed by using three techniques 
to determine the magnitude of lags in the pneumatic angle-of-attack and angle-of- 
sideslip sensing systems over a variety of flight conditions. The techniques were 
applied to stability and control maneuvers and provided comparable results . 

An apparent trend with Mach number showed that the lag in both sensor systems 
increased from approximately 0.15 second at subsonic speed to 0.4 second at Mach 3. 
Because Mach number was closely related to altitude for the available flight data , 
the individual effects of Mach number and altitude on the lag could not be separated 
clearly . However , the results indicated the influence of factors other than simple 
pneumatic lag. 


Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif,, September 9, 1974 
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APPENDIX A 


EQUATIONS FOR TIME-VECTOR ANALYSIS 


The basic longitudinal equation used in time-vector analysis is 

“ = q " V (a n ' 

Expressed in time-vector format , the equation becomes 

l(W 900 + ^) =1 -sItKa. 


where 


fl> d = sin 


The correction for probe displacement from the center of gravity is given by the 
equation 


i „ i a. x 

IflWT 

The preceding equations are solved for the calculated 4> , , which is then used in 

the expression “r q 

- W - ' — Verto* 

The basic lateral-directional equation is 


p = -r + | q> + pa Q + ^ a t 


Expressed in time-vector format , the equation becomes 


>n| 1 1 (*p/r * 90 ° + ♦d) * + f |f |%/r + «o|?|' *p/r + f lrl*.,/r 


The tp/r term can be simplified by using the equations 

IsU-Ufil 


r I w n jr 
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The correction for probe displacement from the center of gravity is given by the 
equation 


<t> 

r i p/r 


d> 


Pj/r 



The lateral-directional equations are solved for the calculated 4> Q 
used in the expression ** 


*«V r >c ' 

u * 555 s 


Period 


j r , which is then 



APPENDIX B 


EQUATIONS USED IN NEWTON-RAPHSON ANALYSIS 


The equations of motion used to determine stability and control derivatives by 
the Newton-Raplison technique of matching flight and calculated data time histories 
are summarized in this appendix . The equations are in body-axis form . 


Longitudinal Derivatives 

The following equations were used to analyze longitudinal short-period 
maneuvers: 


q = M a + M q + M K 6 
n a q o e 

4 o 


q-N 6 6 
^ e o 


v 

a n-T(V + V*V e ) 

8=q 

Position correction for instrument location was unnecessary for a^ , because 

the accelerometer was at the center of gravity; however, the following equation was 
used to correct for the o-probe location: 


a 


a. 

1 


+ 



q 


Lateral-Directional Derivatives 


The analysis of lateral-directional maneuver derivatives was based on the 
following equations: 


P 


I 


f"v 


’xz. N _ 

Y p + V 


+ v + L pP + L 6 8 a + L 6 8 r 
r a r 

•V'V"! 6 a + N 6 8 r 

a r 
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M-r*f?*p. 0 *Yj 1( Y r r ( 0*V ! V Y 6 8 r 
o a r 

v n 

a A = -—(Y p + Y r + Y„0 + Y- 6 + Y e 6 \ 
t gl p F p 0 P 8 fi a ® P v 

9 = p + r0 o 

Position correction for instrument location was not necessary for a { because 

the accelerometer was at the airplane center of gravity; however, the following 
equation was used to correct for the 0-probe location: 






APPENDIX C 


EQUATIONS USED IN THE LEAST-SQUARES 
ERROR MINIMIZATION TECHNIQUE 


To calculate the a and (3 time histories from available accelerations and 
angular rates , the following two equations were used in the least-squares 
technique: 


® = q "V (a n" 1} 


P = -r + a Q p + $ (<p + a t ) 

To obtain traces of a and p , the two parameters were integrated over time for the 
duration of the maneuver at the chosen sample rate . 

Unknown parameter biases were accounted for in the analysis. 



APPENDIX D 


CALCULATION OF CLASSICAL THEORETICAL 
TIME CONSTANT CURVES 


The following representative expression for pneumatic lag in pressure instru- 
mentation is given in reference 1: 

_ 32p (Tubing length) 2 f. + Instrument volume 1 

" (Tubing diameter) 2 yP L (Tubin e len S th > (Tubing cross-section area)J 

Most expressions of this type from other references (refs . 2 and 3 , for example) 
follow this general form and can be summarized as: 



The constant K is some function of y and instrument dimensions and would not 

7 

change for any given installation . For the a system, K = 10.9X10 , and for the 

7 

P system, K = 8.44 X 10 . (For this simplified analysis, only the 640-centimeter 
(252-inch) piece of tubing and the sum of the transducer volumes were used.) The 
pressure, P , is the actual pressure , P 4 g 0 , that was sensed at the a and p ports. 



recorded. However, wind-tunnel tests v"»re made on the complete pitot-static and 
hemispherical head installation to statically calibrate the a and p measurements. 
Figure 15 summarizes these data as the ratio of port pressure to stagnation pres- 

^ 45 ° 

sure , -= — , versus Mach number . The figure also presents some limited flight 
*2 

data obtained with a similar hemispherical head installation on a different airplane 
(ref. 14) . The agreement between the flight data and wind-tunnel data is good. 
Also plotted is the pressure ratio obtained from modified Newtonian theory , which 
can be expressed as 




0 + cos 


2 


0 


where 0 = 45°. 

The overall agreement of the pressure ratios in figure J 5 provides confidence 
in the ability to calculate port pressures for the actual flight conditions of the data 
in this report . 


i 
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• Flight data (ref. 14) 
o Wind tunnel. a = 0 1 ’. , 

Reynolds number = 62.0 x io 6 per m 



Figure 15. Comparison of wind-tunnel-determined pressure 
ratio with theory and flight data. 


Viscosity , p, is a function of temperature only and can be calculated by using 
the following expression from reference 15: 

In SI Lnits , 

2 

H = 145.82 T - \ l0 X 10~ 8 
and in U.S. Customary Units, 

2 

)‘ = 2 - 270 TTT9O XI0 ' 8 

However, the temperature varied throughout the tubing, so an average weighted 
value was necessary. Approximately 85 percent of the tubing was inside the 
YF-12A nose cone, which was cooled to an estimated temperature of 311° K (560° R) . 
The temperature in the nose boom portion of the tubing ranged from the estimated 
311° K (560° R) at the boom support to approximately total temperature at the 
stagnation point of the hemispherical head . 

A theoretically precise determination of the lag was not possible because of 
unknown temperature gradients and the lack of actual port pressures. However, 
some insight into the applicability of the classical theory was obtained by calcu- 
lating port pressures based on wind-tunnel and flight data and using an average 
weighted temperature of 333° K (600° R) . Although an estimate of the temperature 
was used , the results were not adversely affected because lag is not extremely 
sensitive to temperature; for example, an increase of 56° K (100° R) would result 
in only a 12-percent increase in lag. 
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